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ABSTRACT 

Q-i We present a study on the determination of rotation periods (P) of solar-like stars from 

the photometric irregular time-sampling of the ESA Gaia mission, currently scheduled 
for launch in 2013, taking into account its dependence on ecliptic coordinates. We ex- 
amine the case of solar-twins as well as thousands of synthetic time-series of solar-like 
stars rotating faster than the Sun. In the case of solar twins we assume that the Gaia 
unfiltcrcd photometric passband G will mimic the variability of the total solar irradi- 
ance (TSI) as measured by the VIRGO experiment. For stars rotating faster than the 
Sun, light-curves are simulated using synthetic spectra for the quiet atmosphere, the 
spots, and the faculae combined by applying semi-empirical relationships relating the 
level of photospheric magnetic activity to the stellar rotation and the Gaia instrumen- 
tal response. The capabilities of the Deeming, Lomb-Scargle, and Phase Dispersion 
Minimisation methods in recovering the correct rotation periods are tested and com- 
pared. The false alarm probability (FAP) is computed using Monte Carlo simulations 
and compared with analytical formulae. The Gaia scanning law makes the rate of 
correct detection of rotation periods strongly dependent on the ecliptic latitude (/3). 
We find that for P ~ 1 d, the rate of correct detection increases with /3 from 20-30 per 
cent at /3 ~ to a peak of 70 per cent at /? = 45°, then it abruptly falls below 10 per 
cent at /3 > 45°. For P > 5d, the rate of correct detection is quite low and for solar 
twins is only 5 per cent on average. 

Key words: stars: rotation - stars: late-type - methods: data analysis - surveys - 
open clusters and associations: general. 



1 INTRODUCTION 

Main-sequence stars with spectral types later than F5 show 
variability phenomena due to solar-like magnetic activity. 
As in the Sun, magnetic fields generated in the convection 
zone produce active regions (ARs) consisting of cool spots 
and bright faculae. The visibility of the ARs is modulated by 
the rotation of the star and the associated modulation of the 
optical flux induced by their brightness inhomogeneities can 
be used to derive its rotation period. Nevertheless, active 
regions are far from being ideal tracers for stellar rotation 
studies because they are subject to an intrinsic evolution on 
different timescales, ranging in the Sun from a few hours to 
several months in the case of the activity complexes. More 
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precisely, lLanza et al.l (|2004l l studied the typical time-scales 
of the solar variations by analysing time-series of the total 
solar irradiance (TSI) and the spectral irradiance at 402, 
500, and 862 nm (SSIs) as observed by the VIRGO exper- 
iment on board of the SoHO satellite from 1996 to 2002. 
They showed that the intrinsic evolution of sunspots occurs 
on a typical time-scale r s ~ 9 d whereas the evolution of the 
faculae has a time-scale of the order 77 ~ 60 d. The vari- 
ability due to the rotational modulation has a time-scale set 
by the solar rotation period and is of the order of 30 d. On 
time intervals up to 60 d the variation of the solar flux is 
therefore dominated by the intrinsic evolution of ARs and 
the rotational modulation. Variations on a longer time-scale, 
up to 200 d, are attributed to the intrinsic evolution of ARs 
complexes at " active longitudes" i.e. heliographic longitudes 
characterised by the frequent, localised emergence of new 
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magnetic flux (e.g.. iDonahue fc Kei]|[l995l ). Finally, at even 
longer time-scales, the flux variations are mainly due to the 
11-yr sunspot cycle. 

The study of variability phenomena in solar-like stars 
of different ages and spectral types is crucial to understand 
how the stellar rotation and the magnetic activity evolve 
in time and how they are related to each other and to the 
global stellar parameters. The analysis of long-term pho- 
tometric time-series of solar-like stars in the optical pass- 
bands, like those that will be provided as part of Gaia mis- 
sion data-flow, is a powerful tool to investigate the variabil- 
ity induced by magnetic activity. The aim is, in the first 
place, to extract rotation periods from the rotational mod- 
ulation and to correlate them with the magn etic activity 



level i nferred by the light - curve amplitude (e.g..lLanza et al 
200f ; [Lanza et all l2006bl ; iMessina fc Guinanll2002l ; Messina 



2008). The rotational modulation can be used to further con- 



strain the geometrical and physical pr operties of ARs (e.g., 
lLanza et alj|200ll ; IMessina et al.ll2006h . Observations during 
different seasons can further reveal information on the sur- 
face differential rotation, since the ARs latitudes change in 
time because of the magnetic cycle. Finally, long-term flux 
variation c an allow us to estimate magnetic cycles' time- 
scales (e.g.. IMessina fc Guinanll2002l ). 

Gaia is an ESA (European Space Agency) mission that 
will perform a multiepoch survey of the whole sky with the 
aim to chart a three-dimensional ma p of the Milky Way 
l|Perrvman et al.ll200ll ; lMignardll2005l ). The mission, sched- 
uled for launch in 2013, will have a duration of about 5 
years and will supply astrometric, photometric, and spectro- 
photometric measurements for about 1 billion sources down 
to the limit magnitude V « 20. A crude estimate suggests 
that about 20 milli ons of the observed sou rces will be de- 
tected as variables (lEver fc Cuvpersll2000l ). Therefore this 
survey offers a great opportunity to broaden significantly 
the statistics and the characterisation of solar-like variables. 
At the end of the Gaia mission a photometric time-series in 
a broad unfiltered passband (the so called G-band) with a 
time baseline of about 1800 d will be available for each source 
observed by the satellite. The number of per-source obser- 
vations and, consequently, the number of points for each 
time-series, depends on the source coordinates and ranges 
fro m about 40 to 250 with a mean value of 80. 

lEver fc Mignardl l|2005h studied the rate of correct de- 
tection of periods with Gaia as a function of the ecliptic co- 
ordinates for strictly periodic variable stars. They performed 
an extensive analysis by simulating thousands of light curves 
characterised by different periods and signal-to-noise ratios. 
They showed that the rate of detection depends mainly on 
the ecliptic latitude of the simulated sources and that, for a 
signal-to-noise ratio S/N ^ 1.5, it is close to 100 per cent 
for all simulated periods. 

The variability simulated by lEver fc Mignardl (|2005l ) 
consists, however, of a simple sinusoidal signal with a fixed 
frequency so their results do not apply to solar-like vari- 
ability. Indeed, light-curves of solar-like stars are more com- 
plex than a simply sinusoidal signal because the amplitude 
and phase of the flux modulation change in time owing 
to t he intrinsic evolution of ARs and AR complexes (see 
e.g. IMessina fc Guinanl [20021 . IMessina fc Guinarj|2003l ). In 
the present work, we simulate light-curves of solar-like stars 
characterised by different rotation periods and magnitudes, 



and analyse Gaia's capability of correctly detecting rotation 
periods. Moreover, we tested the efficiency of different pe- 
riod search algorithms in order to find out the method that 
gives the best performance for this kind of variability given 
the peculiar Gaia sampling. In Sect. 2 we give a brief de- 
scription of the Gaia instruments and illustrate the main 
features of the Gaia scanning law. In Sect. 3 we describe 
the methods used to simulate Gaia photometric time-series 
for solar twins and for solar-like stars with a rotation period 
ranging from 0.3 d to 180 d. In Sect. 4 we discuss the re- 
sults obtained by running the period search algorithms on 
the simulated time-series. 



2 THE GAIA MISSION 
2.1 The instruments 

Gaia will perform its observations from a controlled 
Lissajous-type orbit around the L2 Lagrange point of the 
Sun and Earth-Moon system. Gaia has two telescopes with 
two associated viewing directions with a field of view of 
0.7° x 0.7° each. The two viewing angles are separated by 
a highly-stable basic angle of 106.5°. The two field of views 
are combined into a single focal plane covered with CCD 
detectors. By measuring the instantaneous image centroids 
from the data sent to ground, Gaia measures the relative 
separations of the thousands of stars simultaneously present 
in the combined fields. 

The astrometric field (AF) in the focal plane is sampled 
by an array of 62 CCDs, each read out in TDI (time-delayed 
integration) mode, synchronised to the scanning motion of 
the satellite. Stars entering the combined field of view first 
pass across dedicated CCDs which act as a sky mapper 
(SM) - each object is detected on board and information 
on its position and brightness is processed in real-time to 
define the windowed region read out by the following CCDs. 
Before stars leave the field of view, spectra are measured 
in three further sets of dedicated CCDs. The BP and RP 
CCDs - BP for Blue Photometer and RP for Red Photome- 
ter - record low-resolution prism spectra covering the wave- 
length intervals 330-680 and 640-1000 nm, respectively. In 
addition to the low-resolution spectro-photometric instru- 
ment, Gaia features a medium-resolution integral-field spec- 
trograph, the so-called Radial Velocity Spectrometer (RVS) 
instrument. 

The astrometric measurements will be performed in the 
G-band, i.e. a broad unfiltered passband that covers the 
wavelength range from about 350 to 1000 nm, with the max- 
imum transmission a t ~ 715 nm and a full width at half 
maximum of 408 nm (|jordi et al.l l2006). The expected pho- 
tometric precision for a single measurement is better than 
2 mmag for targets brighter than G = 15 and about 10 mmag 
for the faintest sources (G ~20). All data will be analysed 
through software packages written in the framework of the 
Gaia DPAC (Data Processing and Analysis Consortium). 



2.2 The Gaia scanning law 

During its 5-year operational lifetime, the Gaia satellite will 
continuously spin around its axis at a constant speed of 
60 arcsec s _1 . Therefore, over a period of 6 hours, the two 
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Figure 1. Number of observation groups (short-term sequences 
of observations closely spaced in time) at different ecliptic coordi- 
nates. 1500 coordinate pairs are considered. Two different regimes 
are clearly visible: the average number of groups is about 50 at 
P > 50° and about 30 for /3 < 40° . 

astrometric fields of view will scan across all objects located 
along the great circle perpendicular to the spin axis. Because 
of the basic angle of 106.5° separating the astrometric fields 
of view on the sky, objects transit the second field of view 
with a delay of 106.5 minutes compared to the first field. The 
spin axis, maintained at a fixed angle £ = 45° with respect 
to the solar direction, is in turn characterised by a slow pre- 
cession (about 4° per day) around the Sun, so that the great 
circle mapped by the telescopes slowly changes in time. As a 
result of the scanning law, the time sampling in the G-band 
is quite irregular and the number of observations (transits), 
for a given source, is a function of the ecliptic latitude. Ba- 
sically, each time-series is made by several groups of ob- 
servations, i.e. short-term sequences of observations closely 
spaced in time. lEver fc Mignardl l|2005h studied the features 
of the Gaia scanning law in details and showed that two 
different regimes can be identified: time-series will consist of 
about 30 groups of observations for sources placed at ecliptic 
latitudes —30° ^ /3 ^ 30°, whereas in the directions above 
P — 40° and below P = —40° the average number of groups 
of transits is 50. Such regimes are clearly visible in Fig. fTJ 
where we plotted the number of observation groups Ng for 
1500 sky directions uniformly distributed in the hemisphere 
delimited by 0° sC A ^ 360°, 0° ^ ,5 ^ 90°. 

The number of observations per group has a more com- 
plicated dependence on the ecliptic coordinates. At latitudes 
P ^ — 45° and P ^ 45° a group consists, on average, 
of only two observations spaced by a six-hours interval cor- 
responding to two consecutive satellite revolutions. At lat- 
itudes —45° ^ /? ==5 45°, the number of observations per 
group ranges between 2 and 60. In Fig. [2] we plot the max- 
imum number of transits per group for each of the ecliptic 
coordinate pairs considered. 



3 SIMULATIONS 

In this section we illustrate the methods used to simulate 
photometric time-series in the Gaia G-band for solar twins 
and for solar-like stars rotating faster than the Sun with 
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Figure 2. Maximum number of observations (transits) per group 
for the same coordinate pairs of Fig.fT] 

different apparent magnitudes. All simulated sources are 
placed at different ecliptic coordinates to study the rate of 
correct period detection as a function of the observing di- 
rection. 



3.1 Simulation of photometric time-series of 
solar-twins 

The total solar irradiance variability closely mimics the 
changes of the solar flux in the optical band (e.g., 
IFrohlich et ai1ll997l . iFligge et al.lll998T ). Therefore, to first 
order, we can assume that the TSI variations, as measured 
by the VIRGO experiment on board of the SoHO satellite, 
can be taken as a proxy for the flux variations as measured 
by Gaia in the integrated G-band for a solar-twin star. The 
TSI time-series used in the present work consists of one mea- 
surement per hour, each measurement having a relative ac- 
curacy of 2.0 x 10 -5 . This time-series is the level 2.0 series 
of TSI made available by the VIRGO team on the official 
web sit^Bof the instrument. The TSI measurements we con- 
sidered span the time interval between 7 February 1996 and 
5 March 2008 and, therefore, overlap with a complete so- 
lar activity-cycle. They show occasional gaps with typical 
durations of a few days, except for the four-months gap in 
mid 1998, when the spacecraft tracking was lost. Details 
on the p rocedure appl i ed to o btain level 2.0 data c an be 
found in lAnklin et all (I1999T ). IFrohlich fc Finsterld (|200lh 
and IFrohlich! (|2003T l. 

For our purposes, we resampled the TSI time-series ac- 
cording to the Gaia scanning law at different values of eclip- 
tic coordinates (A, P) to simulate Gaia time-series. We chose 
770 directions uniformly distributed in the sky region de- 
limited by 0° ^ A < 180° and 0° P 90°. We restricted 
our simulations to ecliptic latitudes 0° ^ P ^ 90° because 
the features of the Gaia sampling are quite symmetric with 
respect to the ecliptic plane. The restriction in longitude 
is justified by the fact that the rate of correct period de- 
tections for strictly periodic variables depends mainly on 

1 http://www.pmodwrc.ch/. The measurements are listed in the 
file virgo.tsi.h_v6.001.0803.dat 
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Figure 3. The TSI time-series as measured by the VIRGO ex- 
periment (black dots). A 5-yr (mission duration) window sampled 
with the Gaia scanning law at (A = 45°, /3 = 45°) is over-plotted 
on the TSI measurements in order to simulate photometric time- 
series acquired by Gaia in the G-band. The window is shifted 
along the TSI time-series to cover different phases of the solar 
cycle (green bullets, red triangles and blue diamonds). 



the ecli ptic latitude and show o nly little dependence on lon- 
gitude (jEver fc Mignardll2005h . For each chosen direction 
(Ao,/3o), we select a time interval of 5 yrs (the duration of 
the Gaia mission) in the TSI time-series and then choose, 
in such an interval, TSI data points according to the Gaia 
scanning law sampling at (Ao,/3o)- The 5-yr window was 
then shifted with steps of 300 days along the entire TSI 
time-series. In this way, for each chosen sky direction, we 
obtained 10 time-series which sample different parts of the 
solar activity cycle with the Gaia time-sampling. In Fig. 
three of the time-series generated at (A = 45°, /3 = 45°) are 
over-plotted on the whole TSI time-series. The simulated 
time-series are characterised by the same number of points 
and time sampling, but cover different intervals of the solar 
variability along the 11-yr cycle. 

The sampling at a given sky position has been gener- 
ated through AGISLab (Astrometry Global Iterative Solu- 
tion Laboratory), a packag e developed by t he Coordination 
Unit 3 of the Gaia DPAC (|Holl et al.ll201ol ). Finally, we de- 
graded the time-series accuracy by adding the instrument 
noise (according to the error model described in lJordi et al.l 
l200fj ) to simulate stellar targets with an average magnitude 
G=10. 

3.2 Simulation of photometric time-series for 

solar-like stars with different rotation periods 

3.2.1 The model 

We used the approach described bv lLanza et al.l l|2006af ) to 
generate synthetic light curves of solar-like stars. It assumes 
that the stellar variability can be modelled by means of three 
discrete active regions plus a un iform background compo- 
nent. This is based on the work of [Lanza et all <|2003l , where 
the TSI time-series is fitted with this simple model with 
a residua l standard deviations o f 30-35 part s-per-million 
(see also iBonomo fc Lanzal 2008). Following |l anza et al.l 



l|2006al ). we simulate the flux variations at a given wave- 
length A and time t as: 

AF(\,t) = F(X,t) - F U (X) = 

= ^2 Li k A k I u (\, Li k ) [c S p t(A) + Qc fac (A,/i fc )] 

k-H k >0 

+AF back (A,i), 

F U (X) is the unperturbed flux (i.e. the monochromatic 
flux coming from the star in the absence of magnetic ac- 
tivity) and F(\,t) the perturbed flux at wavelength A and 
time t; A k is the area, at time t, of the cool spots asso- 
ciated with the k-ih AR (k —1,2 or 3). The contribution 
of surface brightness inhomogeneities to the stellar flux is 
proportional to their projected area over the stellar disk, so 
A k is multiplied by the factor fj, k = cos tf)k , where tf)k is the 
angle between the normal to the k-ih AR, assumed to be 
point-like, and the line of sight; J u (A,/i) is the specific in- 
tensity of the unperturbed photosphere, c spot the contrast 
of cool spots, Cf ac the facular contrast, Q the ratio of the 
facular area to the spotted area in each active region (as- 
sumed to be constant) and AFb ac k(A,i) the flux variation 
due to the uniformly distributed background. The contrast 
factors c S pot and Cf ac depend on AT spot and AT fac respec- 
tively, where AT spo t is the average temperature difference 
between the unperturbed photosphere and the spotted ar- 
eas while ATf ac is the temperature difference between fac- 
ulae and the unperturbed pho tosphere (for detai ls on the 
evaluation of c apo t and Cf ac , see lLanza et al.ll2006af ). 



3.2.2 Simulation of synthetic spectral time-series for 
solar-like stars 

Time-series of synthetic spectra for solar-like stars were sim- 
ulated ac cording to equation CO . The method, described in 
details in lLanza et a l. (2006a), provides spectra F(\,t) of a 
solar-like inhomogeneous photosphere by combining specific 
intensities 7(A, n) computed assuming plane-parallel model 
atmospheres with effective temperature (Teff) equal to that 
of the unperturbed photosphere as well as with T e s differ- 
ing from the unperturbed photosphere by AT apo t and ATf ac 
in the spots and faculae, respectively. The specific inten- 
sities used in this work were computed using of MARCS 
|Gustafsson et all 120081 ) for logg = 4.5 and solar chem- 
ical composition. The parameters (areas and coordinates) 
of the three active regions and the background component 
AFb a ck(A,t) at a given time t are computed by interpolat- 
ing in time the values supplie d from the best-fits of the TSI 
time-series (|Lanza et al.l l2003). The dependence of the rota- 
tional modulation amplitude on the rotation period P and 
T e ff is taken into account by scaling t he spot areas and th e 
background component as described in lLanza et al.l(|2006al ). 

As illustrative examples, in Fig. [4]we compare the sim- 
ulated spectrum of a magnetically active star with the cor- 
responding unperturbed spectrum and in Fig. [5] we show 
the effects on the light-curve of varying Q and P. Note that 
in stars whose variability is "faculae dominated", the rota- 
tional modulation is less evident because faculae produce an 
increase of the optical stellar flux that counteracts the effect 
of cool spots (Fig. O top panel). In stars rotating faster the 
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Figure 4. Comparison between spectra of an unperturbed 
(black) and a magnetically active (red line) photospheres. 

variability amplitude is larger (Fig. [5] bottom panel) because 
the magnetic activity is higher. 
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Figure 5. Top panel: simulated light curves for two stars having 
the same effective temperature T c g and the same rotation period 
but a different ratio Q between facular and spotted areas. Bottom 
panel: two synthetic light curves for stars with the same effective 
temperature T c g and the same AT spo t and Q values are plotted. 



3.2.3 Simulation of G-band time-series 

The procedure we follow to generate a G-band time-series 
for a star placed at ecliptic coordinates (Ao,/3o) consists of 
the following steps: 

• we run the AGISLab code and compute the transit 
times at (Ao,/3o) according to the Gaia scanning law; 

• for each transit time we compute the spectrum F(X, i) 
as given by equation JT} ; 

• we convolve the spectrum F(\, i) with the transmit- 
tance profile of the Gaia G band to obtain the stellar flux 
in the G-band; 

• we add noise to the light c urve according to the error 
model described in ljordi et alj (|2006t ). 

We simulated photometric time-series in the G-band for 
about 1500 ecliptic coordinate pairs uniformly distributed 
in the region 0° < A < 360° and 0° < f3 < 90° . 
The parameters used for our simulations are: 

• Rotation period P: 0.3, 0.5, 1.0, 2.0, 4.0, 5.0, 7.5, 10.0, 
20.0, 40.0, 60.0, 90, 180 d; 

• Effective temperature T e a= 5000 K; 

• Q =1; 

• AT spot = -800 K; 

• AT fac = 110 K; 

• Apparent magnitudes G= 10, 15, 19. 

In total, we simulated about 58000 photometric time- 
series corresponding to ~ 4.6 x 10 6 transits. 



4 RATE OF CORRECT PERIOD 
DETECTIONS WITH GAIA 

In the present section we show the results obtained by run- 
ning three period search algorithms on the simulated time- 
series. Th e algor ithms we tested are: the Deeming method 
dDeeminsl 19751), the Lomb-Scargle periodogram l|Lombl 
1 19761 : Scargle 19821) and the Phase Disp e rsion Minimizatio n 
method (PDM method: |jurkevichlll97ll : IStellingwerJll978l ). 



We preliminary discuss a method to estimate the FAP (False 
Alarm Probability) for the periods detected by the Lomb- 
Scargle algorithm which is the method providing the highest 
rate of correct detections. 



4.1 FAP estimation 

The FAP associated with a given peak with a power level z in 
the Lomb- Scargle periodogram is the probability that such a 
z value is due to pure Gaussian noise, which is the main error 
source for stars brighter than G — 16 mag in the Gaia pho- 
tometry. In the general case M onte Carlo simulation s can be 
used to estimate the FAP (e.g. jFrescura et al.ll2008l ). These, 
however, are not suitable for the analysis of a very large 
dataset because of the computational load implied. The 
Lomb-Scargle algorithm, however, al lows an analytical esti- 
mate of the FAP which, according to ISchwarzenberg-Czernvl 
l|l998l ). is given by: 



Pr\Z n 



> z] = 1 - 



N/2 J 



N/2' 



(2) 



for a given frequency v with a power z in the normalised 
Lomb-Scargle periodogram. Equation ((2} gives the probabil- 
ity that the maximum peak power Z max exceeds a threshold 
z for a data set {xi} consisting of pure Gaussian noise; N is 
the number of data points used to compute the periodogram, 
cr^ the variance of the data set {xi} an d M the number of 
independent frequencies, i.e. the number of frequencies at 
which the periodogram powers are independent variables. If 

1, equation becomes 



one normalises such that a 2 ^ 



Pr[Zma X > Z] = 1- (l-(l-2 2 /A') iV/2 ) J 



(3) 



Note that, as indicated by ISchwarzenberg-Czernvl 

( 1998 ), only in the limit N — > oo, equation ([2]) reduces to 
the I Home fc Baliunasl (|l986l ) expression: 



Pr[Z n 



> z] 



(1 



(4) 



6 E. Distefano et al. 



which is the distribution often adopted by others authors in 
this contest. 

The difficulty in applying Eq. (J2j and Eq. (j4]) arises from 
the unknown parameter M. In the case of a set of data 
with an even time sampling, the independent frequencies 
are given by 

V K = ji, (5) 

where T is the time interval spanned by the data and 
k = 1, 2, ...N/2. In this case M = N/2. In the case of a set 
of unevenly sampled data, M cannot be derived analytically 
and can be inferred only by fitting empirical cumulative dis- 
tribution functions (CDFs) of the peak powers z generated 
by means of Monte Carlo s imulations fsee lHorne fc Baliunasl 
1986, lFrescura et al.ll2008h . One could expect that M satis- 
fies the condition N/2 ^ M ^ Nf, where Nf is the number 
of inspected frequencies. One possible choice is to adopt the 
value 

M = — (6) 
r 

where r = (l/T)/8v and 5v is the frequency step used to 
sample the periodogram |Frescura et alj [20081 ). The factor 
r is introduced to take into account that a periodogram is 
generally oversampled (i.e. sampled with a frequency step 8u 
smaller than the minimum frequency interval Ai/ m i n — 1/T 
that can be resolved by the dataset). In our tests, we com- 
puted the periodogram in the frequencies domain (0.01 - 5 
d" 1 ), with a frequency step 5v = 0.0001 d" 1 and, there- 
fore, Nf = 49990, Au min = 1/1800 d" 1 = 0.00055, r = 5 
and M = N f /r = 9980. We used Monte Carlo simulations 
to construct an empirical FAP for 8 different sky positions 
characterised by a different sampling and number of obser- 
vations and made a comparison between such empirical dis- 
tributions and equation ((3J computed adopting M — 9980. 

The approach we followed to generate the empirical 
FAP at a given pos ition (Ao,/?o) ' s the same described in 
iFrescura et ail (|2008r i. i.e.: 

• we generated 10 4 time-series of pure white Gaussian 
noise with the sampling given by the Gaia scanning law at 
(Ao,A>); 

• we run the Lomb-Scargle algorithm on each time-series 
by using the same grid of frequencies applied to period 
search in our tests; 

• for each of the 10 4 periodograms we selected the high- 
est power Zmax and used these values to construct their 
cumulative distribution functions (CDFs). 

A CDF constructed in such a way can be regarded as 
an empirical representation of the probability function 
Pi[Zmax £ z ]y i- e - the probability that pure noise alone 
could generate a peak greater than or equal to a given 
threshold value z in the Lomb-Scargle periodogram. In Fig. 
[U]we show the CDFs constructed for five different sky posi- 
tions. The CDFs depend on the number of observations, i.e. 
the lower the number of observations the higher the proba- 
bility of finding high z peaks in a periodogram generated by 
a set of data consisting of pure Gaussian noise. In the same 
figure we also plotted the CDFs constructed for two differ- 
ent sky positions, (A = 0°,/3 = 75°) and (A = 160°, ,5 = 0°), 
in which the number of observations is the same (TV = 80) 
but the sampling is different, the first being more uniformly 



distributed, the second concentrated in 12 rather isolated 
groups. The two CDFs are fairly similar, which indicates 
that the CDFs are much more sensitive to the number of 
observations than to the observations' distribution in time. 

In Fig. we compare the Monte Carlo CDFs with the 
probabilities given by equation ^ and the number of inde- 
pendent frequencies estimated by equation©, for two dif- 
ferent sky directions. The comparison shows that the FAP 
obtained using equations © and (O underestimates the true 
FAP. In Fig. [7] we also plot the results obtained by fitting 
both equations ([3| and ((4]) to the empirical CDFs using a 
non-linear least-squares (NLLS) Marquardt-Levenberg algo- 
rithm, varying M as a free parameter. These best-fits indi- 
cate that the number of independent frequencies estimated 
using equation © is too low and that a higher value should 
be adopted. The shape of the CDFs is not accurately re- 
produced either by equation © or by equation Q, since 
the former tends to underestimate the FAP while the latter 
tends to overestimate it. 

In conclusion, an accurate FAP estimate would require 
the construction of CDFs for a convenient grid of sky co- 
ordinates. However as a first approximation, we shall use 
equation ((3]) adopting a suitable number of independent fre- 
quencies M. To do that, we fit equation© to the 8 em- 
pirical CDFs (corresponding to the selected 8 different sky 
positions) in the region Pr < 0.1 (corresponding to a signif- 
icance level > 90 per cent), then we take the average of the 
resulting M values, which gives Mad =< M >= 27500, our 
adopted number of independent frequencies. In Fig. [8] we 
display the comparison between empirical CDFs and equa- 
tion ([3]) computed by adopting M = Mad- The empirical 
CDFs are not accurately reproduced in the region where 
FAP > 0.1, but are reproduced with an accuracy better 
than ~ 5 per cent in the region where FAP < 0.1 which is 
the region of interest for all practical purposes. 

4.2 Search for rotation period in solar-twins 

In Fig. O we show the distributions of periods detected 
by using the Lomb-Scargle algorithm on solar-twin simu- 
lated time-series with a FAP < 10~ 2 , 10" 4 , and 10~ 5 . For 
FAP < 10 -2 , the true period (i.e. the solar synodic period 
Pq = 27.27 d) is recovered, with a relative error less than 
10 per cent, only in a small fraction (about 5 per cent) of 
the 7700 simulated time-series, while in the 60 per cent of 
cases periods shorter than about 20 d are obtained. By de- 
creasing the FAP threshold from 10~ 2 to 10" 4 , we obtain 
a minor decrease in the percentage of the recovered true 
period and a large decrease of the percentage of spurious 
periods below 20 d. These latter decrease below 15 per cent 
for FAP < 10" 5 , at the expenses of a diminished true period 
recovery at 3 per cent. 

The results shown in Fig. [9] are not surprising. In solar- 
like stars, the variability due to rotational modulation can 
be masked by the growth and decay of ARs and AR com- 
plexes. Such variability phenomena are responsible for the 
other periods detected by the period search algorithms as 
shown in the distribution. As we pointed out in the intro- 
duction, the time-scale of intrinsic evolution of a single AR 
in the Sun is comparable with the solar rotation period, but 
AR complexes at activ e longitudes have a lifetime of about 
t ~ 200 - 250 d l|Lanza et al.l l2003h . In other words, 
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Figure 6. Left panel: empirical CDFs constructed for five different sky positions (A,/3). Right panel: enlargement of the critical region 
Pr[Z < z] < 0.1. Note that the CDFs relative to two different sky positions, (A = 0° , ft = 75°) and (A = 160°, (3 = 0°), with the same 
number of observations but differently distributed appear to be very similar. 




Figure 7. A comparison between empirical CDFs (black dotted lines) and the analytical expression in equation J3]) (in purple) computed 
with the number of independent frequencies evaluated using equation foil. Also plotted the results obtained by fitting the CDFs with 
equations J3j) (in green) and iflj (in red) by varying M as a free parameter. 
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Figure 8. Left panel: a comparison between empirical CDFs (black dotted lines) and the analytical expression J3jl (in purple) computed 
by using M = M a j = 27500. Right panel: enlargement of the critical region Pr[Z < z] < 0.1. 
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Figure 9. Distributions of the periods as detected by Lomb- 
Scargle algorithm on simulated G-band time-series for a solar- 
twin star when time-series of 5-year extension are analysed. 



even though a single spot group has a short lifetime (one 
or two weeks), AR complexes may persist for several solar 
rotations because new spot groups form as previous ones 
fade away. Therefore, the solar flux variations due to the 
evolution of AR complexes mask the variation s due to the 
rotati on on time-scales longer than 200 — 250 d. lLanza et al.l 
(5553) showed that the rotation period of the Sun can be 
recovered only if long time-series are divided into segments 
with a time extension ~ 150 d to limit the effect of AR 
complexes evolution. Moreover, the period was correctly de- 
tected only on segments falling close to the minimum of 
the solar cycle. In fact, during such a phase, the solar pho- 
tosphere is dominated by faculae, which evolve on a time- 
scale Tf ac ~ 60 — 80 d ~ 2 — 3 Pq. Therefore, close to 
the minimum of the 11-yr solar cycle, the signal induced by 
rotational modulation is more coherent than in the inter- 
mediate phase and the maximum, when the level of solar 
activity is higher and sunspots dominate the TSI variations. 
The sunspot evolution on a time-scale of 10-15 d, on the 
other hand, makes it difficult to detect the rotational modu- 
lation on a time-scale of 27 d. Variability with a periodicity 
of the order of half the solar synodic rotational period may 
be due t o the (still deba ted) presence of active longitudes 
(see, e.g. [Balthasar 2007, and references therein), but the 
frequency of periods around this value in Fig. [9] even for 
FAP < 10 _J , is too high to be consistent with the spectral 
power (relative to the synodic period) obtained by specific 
studies not limited by the Gaia sampling. Therefore, the rel- 
atively high frequency of periods below 20 d are likely due 
to the combination of the lack of coherence of the signal and 
the peculiarities of the Gaia sampling. Note also the pres- 
ence of a persistent peak at periods around 58 d, which is 
not associated to any known variability time-scale. Further 
detailed analysis on these issues is required but deferred to 
a future work. 

The Deeming and PDM method give approximately the 
same rate of correct detections and are not discussed in de- 
tails here. 
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Figure 10. A 160 d sub-series extracted at (A = 45°, = 45°). 
In spite of the high number of observations, the flux modulation is 
poorly sampled because most of the points are distributed along 
an interval of only 7 d. 



4-2.1 Time-series segmentation for solar-twins 

Keeping in mind the results of lLanza et alj (|2004l ) . we tried 
to increase the rate of correct period detection by splitting 
each of the 7700 simulated time-series into sub-series (seg- 
ments) with a shorter extension r. 

The choice of r is quite critical because it should be 
of the same order of the time-scale of AR complexes evo- 
lution but, at the same time, long enough to include a suf- 
ficient number of data points. In our case, to detect the 
rotational modulation of a slowly rotating star like the Sun, 
the length of the sub-series should be selected in order to 
maximise the number of groups Ng per sub-series rather 
than the number of single observations N. In Fig. [TU] we 
show a sub-series, extracted from one of the simulated time- 
series at (A = 45° , /3 = 45°), with a basel ine r = 160 d 
(the same value used by lLanza et al. I bin ll i and consisting 
of a rather large number of single observations (N — 55) 
but mostly concentrated on a single group spanning a 7-d 
interval. Clearly, such a kind of sampling is inadequate to 
trace the solar flux modulation with a period Pq = 27.27 
d. 

The typical interval from two successive groups is about 
45 d (|Ever fc Migna rd 2005), thus we decided to use a base- 
line t = 450 d in order to generate sub-series including, 
on average, 10 groups of observations each. We split each 
of the 7700 simulated time-series into 450 d segments and 
run the period search algorithms on all sub-series satisfying 
the condition Ng 10. We processed about 110000 sub- 
series and plotted the distributions of detected periods in 
Fig. 1111 Unfortunately, in our case, time-series segmenta- 
tion has the effect to reduce the rate of correct detections. 
We also tried to process th e sub-series with the string-length 
method l|Dworetskv II 19831 ). that is particularly suitable for 
dealing with time-series with a low number of points. Such 
a method is based on the minimisation of a quantity, called 
string-length, that is the sum of the lengths of the segments 
joining two successive points in the phase diagram. This 
method is slightly more effective than the others but the 
rate of correct detections is always of the order of 5 per 
cent. 
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Figure 11. Distribution of detected periods after applying time- 
series segmentation 
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Figure 12. Rate of correct detections with the string-length 
method: the top-panel refers to the case in which the solar period 
is detected as the first minimum of the periodogram whereas the 
bottom panel refers to the case in which the solar period coincides 
with one of the three periods with the highest probability. 



Finally, we analysed how the percentage of correct de- 
tections depends on the ecliptic coordinates of the simulated 
source. In Fig. 1121 we show the results obtained by running 
the string-length method on sub-series. For each of the sim- 
ulated sky directions, the percentage of sub-series in which 
the rotation period of the Sun has been recovered with a rel- 
ative error $C 5 per cent is plotted. The picture shows that 
the rate of correct detections (that has an average value ~ 5 
per cent) is not uniformly distributed and peaks at ecliptic 
latitudes close to /? = 45° where it reaches ~ 15 — 20 
per cent. The detection rate obtained by running the other 
period search methods on sub-series have a trend similar to 
that shown in Fig. 1121 and are not shown here. 

The main results of our analysis on solar twins can be 
summarised as follows: 

• the rate of correct period detection for solar-twins is 
generally very low and has an average value ~ 5 per cent; 

• segmentation of sub-series does not significantly in- 
crease the rate of correct detections; 



• at ecliptic latitudes close to j3 = 45° the rate of correct 
detections rises up to 15-20 per cent. 



4.3 Search for rotation period in solar-like stars 

In this section, we investigate the Gaia's capability to detect 
rotational modulation in solar-like stars with different rota- 
tion periods. For doing that, we simulated photometric time- 
series in the Gaia G-band for solar-like stars placed at differ- 
ent ecliptic latitudes as described in Sect. l3.2l with a rotation 
period ranging from 0.3 to 180 d an d processed then with 
the Deeming method (IDeemingl 1975 L the Lomb-Scargle pe- 



iodogram (|Lombll 19761 ; Scargle|[l982 ) and t he Phase Disper- 



sion Minimization method (PDM method: |jurkevicblll97ll : 
IStellingwerj|l978l . with 5 bins in the phase diagram). 

The histograms shown in Fig. [13] visualize, for each pe- 
riod search method, the percentage of cases in which the 
rotation period of the simulated stars is correctly detected. 
Histograms in the left panel refer to the cases in which the 
rotation period is measured with a relative error ^ 5 per 
cent whereas those in the right panel refer to the cases in 
which the rotation period is derived with a relative error 
^ 10 per cent. This percentage decreases further with de- 
creasing relative errors limit below 5 per cent (not shown 
in figure), which indicates that, in our case, noise causes 
a rather wide spread of measured period around the true 
value. As discussed in Sect. 4.3.1, the active region evolu- 
tion also contributes to the period spread around the true 
value (see Fig. I15|l . 

The main results of our analysis are: 

• The rate of correct detections decreases with increasing 
rotation period and ranges from about 60 per cent for the 
fastest rotating stars (P = 0.3 d) to about 5 per cent for 
stars with a rotation period P ^ 20 d; 

• The Deeming and Lomb-Scargle methods are more ef- 
fective than the PDM-Jurkevich algorithm. The rotation pe- 
riod P = 0.3 d, for example, is recovered in about 60 per 
cent of the simulated sky directions by running either the 
Lomb-Scargle or the Deeming algorithm, whereas it is re- 
covered only in « 40 per cent of the cases with the PDM- 
Jurkevich method; 

• For stars with a rotation period P ^ 5 d, the rate 
of correct detection does not depend on the apparent mag- 
nitude. In fact the amplitude of the light modulation is an 
increasing function of the angular velo city and, for P ^ 5 d , 
it ranges between 0.0 4 — 0.3 mag (e.g. iMessina et alj|2010al : 
M cssma et al.ll2010bl ). Therefore the rotational modulation 
can be detected also at G ~ 19 where the photometric 
error is ~ 20 mmag. 



4-3.1 Time-series segmentation for fast-rotators 

Our knowledge on time-scales of AR growth and decay 
(ARGD) in solar-like stars in general is still quite limited. 
iDonahue et al.l |l997al ). from the analysis of 35 lower-main- 
sequence stars observed at Mount Wilson Observatory, es- 
timated the time-scale of active region evolution is approxi- 
mately 50 days, while the lifetime of active region complexes 
is on the order of one year. Most of the stars in their sample 
have ARGD time-scales significant ly longer than the rota- 
tion periods. IDonahue et al.l l|l997bl ) extended such analyses 
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Figure 13. Rate of correct period detection on simulated light curves: the percentage of detections is a decreasing function of the 
rotation period; different colours refer to different period-search methods as labelled. The rate of detections at G=10, 15 and 19 mag 
are shown in the top, medium and bottom panels respectively. The left panels refer to the percentage of cases in which the period is 
correctly detected with a relative error ^ 5 per cent whereas the histograms on the right panel indicate the percentage of cases in which 
rotation period is measured with a relative error ^ 10 per cent 



to 100 stars, showing that the variability of older stars (with 
lower activity and slower rotation) tend to be AR-evolution 
dominated, while younger stars (more active and fast ro- 
tating) have AR evolution time-scales significantly longer 
than the rotation period. The age at which this transition 
ta kes place also depend s on the stellar colour (see Fig. 5 
of iDonahue et aLlll997bl ) and it is around the Sun's age or 
older for stars with (B — V) similar to that of the Sun. Subse- 
quently, iMessina fc Guinaa |2003l ) analysed six young solar 
analogues and found that for four of them the time-scale 
at which the evolution of active regions begins to affect the 
variance of the observed flux ranges from ~5 months (~30 
rotation periods) to ~16 months (~180 rotation periods). 
Two of them were found AR evolution-dominated, in which 
case the active region evolution is not sufficiently distinct 
from the rotation period. 

For the analy sis of rotation periods in sta rs of young 
loose associations, M essina et al. | (|2010h| , |201ll ) adopted a 
time-series segmentation of less than 60 days and used the 
comparison of the periods obtained in different segments to 
estimate the confidence on the resulting period. The num- 
ber of segments in which the period was found with a con- 
fidence level higher than 99 per cent and/or the number 
of segments in which the phased light-curve shows average 
residuals lower than its amplitude was used to estimate the 
robustness of the period determination. Their approach has 
allowed to establish with very high confidence the rotation 
period for ~60 per cent of the entire catalogue (290 stars in 
total). 

It is therefore expected that segmenting the whole Gaia 
time-series would be more effective for young, rapidly ro- 
tating stars, in which the ARGD time-scales are in general 
significantly longer than the rotation period. 

In Fig. [TJ we plot a simulated light curve of a solar- 
like star with a rotation period P — 2 d and superimposed 
on it the points as sampled by the Gaia scanning law at 
(A = 45°, /3 = 45°). The sequence of 49 points refers to the 
same sampling of the observations shown in Fig. 1101 The 
comparison of the two figures shows how the Gaia scanning 
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Figure 14. A simulated light curve for a solar-like star with 
a rotation period P = 2 d (solid line); points corresponding 
to the longest observational epoch at (A = 45°, = 45°) arc 
superimposed on it (red bullets). 



law favours (at least at low ecliptic latitudes) the detection 
of the rotational modulation in fast rotating stars. Note that 
in Fig. [10] the sequence of 49 points covers only one fourth 
of the Sun rotation period whereas in Fig. [14] it traces three 
consecutive stellar rotations without the daily gaps typical 
of ground-based observations. 

It is therefore evident that running the period search 
algorithms on a single group of observations like that shown 
in Fig. [T?] rather than on the whole time-series would be 
more effective because the signal is more coherent along such 
a short-term sequence than in the whole 5-yr time-series. 

We therefore performed a test to verify if, for rapidly 
rotating stars, the rate of correct detections can be increased 
by processing such short sub-series. We selected simulated 
time-series for stars with rotation periods P = 0.3 d and 
placed them at ecliptic coordinates where there is at least 
one group of observations with N ^ 10: in such a way 
we obtained a sample of 377 time-series. Then, we extracted 
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Figure 15. Red histogram: distribution of periods detected by 
running the Lomb-Scargle algorithm on the whole time-series. 
Green histograms: period distribuiton obtained by running the 
Lomb-Scargle algorithm on sub-series. Both distributions are 
centered on the correct period but the distribution in green is 
sharper, indicating a more accurate measure of the rotation pe- 
riod. 



the groups of observations satisfying the condition N > 10 
from each time-series and ran the Lomb-Scargle algorithm 
on such sub-series. In Fig. [I5]we plot the distribution of the 
periods detected by running the Lomb-Scargle algorithm on 
the whole time-series (histogram in red) as well as the distri- 
bution obtained by running the algorithm on the sub-series 
(in green). Both distributions are centred on the correct pe- 
riod but the distribution plotted in green is sharper than 
the other one. Therefore, running the Lomb-Scargle algo- 
rithm on the short sub-series gives a more accurate period 
estimation. 

It is interesting to notice that about 200 of the 377 
analysed time-series have two or more groups of observations 
satisfying the requirement N ^ 10. For fast rotators placed 
at these sky positions, therefore, 2-5 measures of rotation 
period separated by at least 6 months in time are available. If 
these stars have a solar-like magnetic cycle and the latitudes 
of their ARs change with time, the measure of the rotation 
period at different epochs could even allow us to de tect a 
possible differential rotation fcf. iBaliunas et al] [l985V 

The longest rotation period that we can detect with 
confidence in the Gaia case along a given sky direction by 
processing short-term sub-series can be roughly estimated as 
follows. If we assume that a period is detectable if at least 
one and half consecutive rotations are sampled, then, along 
a given direction (Ao,/?o), the upper limit for a detectable 
period is P max ~ I / 1.5 d where L is the duration in 
days of the longest group of observations at (Ao, j3o). In 
Fig. [16] we plot P m ax as a function of ecliptic coordinates. 
The figure shows that very fast rotators (P ^ Id) can 
be characterised in about 25 per cent of the simulated sky 
directions. At ecliptic latitudes close to /3 = 45°, rotation 
periods up to ~ 5 d can be recovered. 
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Figure 16. The longest detectable rotation period versus the 
simulated sky directions. 
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Figure 17. Rate of correct period detections as a function of the 
ecliptic latitude /3 and the simulated frequency. 



4-3.2 Dependence on the ecliptic latitude 

In Fig. 1171 we plotted the rate of correct detection as a func- 
tion of the ecliptic latitude for some of the simulated periods. 
The periods used corre s pond to the frequencies used in Fig. 
10 of lEver fc Mignardl (|2005r i. thus a comparison between 
our wo rk anf theirs can be don e. 

In lEver fc Mignardl (|2005h . the rate of correct detection 
of strictly periodic signals is very high for all the simulated 
frequencies and almost independent of the ecliptic latitude: 
for a signal-to-noise ratio S/N = 1.3 it is close to 100 per 
cent for all the frequencies and for S/N = 0.9 it is greater 
than 40 per cent. On the other hand, in the case of solar-like 
variability, the rate of correct detection is always lower than 
60 per cent and is very sensitive to the ecliptic latitude. For 
example, Fig. [T7] shows that for stars with a rotation period 
P — 0.3 d, the rate increases with ecliptic latitude in the 
range 0° ^ f3 ^ 45° ; it is about 40 per cent at low ecliptic 
latitudes and reaches about 100 per cent at = 45° . Finally, 
it abruptly falls to about 20 per cent at ecliptic latitudes 
/3 ^ 50°. The detection rate is also very sensitive to the 
simulated periods and for slow rotators (P ^ 10 d) it is al- 
ways l ower than 20 per cent. The results o f lEver fc Mignardl 
l|2005h do not apply to solar-like variability because the am- 
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plitude and the phase of the flux modulation change with 
time with the evolution of the active regions and the level of 
magnetic activity. Such a behaviour is clearly shown by the 
simulated light curves plotted in Fig. [5] where the depth of 
the minima changes with time. These light curves are more 
complex than a simple sinusoidal signal and two points hav- 
ing the same phase may have different amplitudes if their 
separation in time is large enough. Ideally, the detection of 
the rotation period for these stars would require a group of 
observations distributed as in Fig. 1141 Such groups, charac- 
terised by at least 10 observations distributed in one or a 
few days, are concentrated at —45° ^ /3 ^ 45° (cfr. Fig. 
[2]), and this explains why the rate of correct detections, for 
fast rotators, increases in this range of latitudes. 



5 THE CONTRIBUTION OF GAIA TO THE 
STUDY OF STELLAR ROTATION IN OPEN 
CLUSTERS 

In the last years, several projects have been carried out to 
monitor late-type stars in open clusters of different ages with 
the aim of measuring their rotation periods and investigating 
the dependence of the stellar rotation on age, spectral type 
and magnetic activity. Such projects, e.g. the RACE-OC 
project (Rotation an d Activity Evolution in Open Cluster) 
dMessina et al.| [2008) or the Monitor project (Hodgkin ct al. 
120061 ). are important to constrain theoretical models of the 
rotational evolution of stars. Although these studies im- 
proved our knowledge of the rotational evolution, ground- 
based monitoring of open clusters suffers from several limita- 
tions. The first is the uncertainty of the stellar ages. The pri- 
mary method to estimate the age of an open cluster relies on 
comparing theoretical isochrones with the colour-magnitude 
diagram (CMD) of the cluster. This method is very sensitive 
to the uncertainty in the cluster distance and, moreover, it 
is affected by the contamination of the CMD by field stars. 
Other limitations are related to the measure of the rotation 
periods, with a bias towards slow rotators (P ^ 10 d), and 
the aliasing due to the day-night duty cycle of ground-based 
observations. 

Gaia will allow us to improve the estimate of open clus- 
ter ages mainly because it will supply very accurate mea- 
surements of the stellar distances. Parallaxes measured by 
Gaia will have an uncertainty of only 0.1 — 0.3 per cent 
for about 10 of the observed sources and an uncertainty 
^ 1 — 3 per cent for about 10 7 sources. This implies, for 
instance, that at the end of the Gaia mission individual par- 
allaxes will be determined with an uncertainty ^ 0.1 — 0.2 
per cent for all stars belonging to the Pleiades, while its av- 
erage parallax for the Pleiades is presently known with an 
uncertainty of about 8 per cent. The availability of paral- 
laxes and proper motions (and also of radial velocities for 
stars down to V ~ 15) will allow us to assess the member- 
ship and obtain uncontaminated CMDs for all the observed 
open clusters, so the expected uncertainty in their absolute 
ages will become 5 per cent. 

Together with distances and ages, Gaia will supply ac- 
curate measurements of the rotation periods for fast ro- 
tating stars belonging to the open clusters in the range 
—45° ^ ft ^ 45° (i.e. for about 50 per cent of the known 
clusters), plus less accurate measurement for longer periods. 
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Figure 18. The longest detecta ble rotation period is plotted for 
each of the open cluster listed in lDias et al.l l|2002l 'l. 

In Fig. 1181 we plotted the ecliptic coordinates of the open 
clusters as listed in the c atalogue "Optica lly visible open 
clusters and Candidates" ()Dias et al.| [2002) using a colour 
scale to depict the longest rotation period detectable P max 
as a function of the cluster coordinates. Note that most of 
the clusters listed in the catalogue have angular dimensions 
smaller than those of the Gaia astrometric fields of view and, 
therefore, the pattern of the scanning law will not change 
significantly for the stars belonging to the same cluster. Al- 
though the information on rotation periods in open clusters 
will not be complete, since accurate periods will be avail- 
able with confidence only for the faster rotators, the num- 
ber of period measurements will be unprecedented. One can 
therefore envisage supplementing such an information with 
ground-based monitoring programs, like the ones mentioned 
above, to build up complete samples of rotation periods for 
a large number of open clusters. 



6 CONCLUSIONS 

In the present paper we study the capability of detecting 
the rotation periods of solar-like stars through the analysis 
of simulated photometric time-series that mimic those of the 
upcoming Gaia mission. 

We test different period search methods and show that 
the Lomb-Scargle algor ithm is the most efficien t. 

We show that the ISchwarzenberg-Czernvl j 19981 ) false 
alarm probability (FAP) analytical formulae do not repro- 
duce accurately the empirical CDFs constructed by means 
of Monte Carlo simulation unless an ad hoc evaluation of the 
number of independent frequencies M is made. We obtain a 
fit to the empirical CDFs using the ISchwarzenberg-Czernvl 
( 1998) formula using M as the fitting parameter with an 
accuracy better than 5 per cent. 

We show that the rate of correct period detections 
strongly depends on the ecliptic latitude of the observed 
source because the time sampling generated by the Gaia 
scanning law exhibits different patterns at different ecliptic 
latitudes. A high percentage of detections are expected in 
the latitude range —45° ^ /3 ^ 45°. The highest rate is 
expected at f3 = ± 45° where the Gaia scanning law gives 
the highest number of transits. 



Solar-like rotation periods with Gaia 13 



The detection rate depends also on the stellar rotation 
period and, in particular, the Gaia scanning law is mostly 
suited to detect periods of rapidly rotating stars (P ^ 5 d) . 
On the other hand, the rate of correct detections for solar- 
twins is, on average, ~ 5 per cent with a maximum of 15 — 20 
per cent at ecliptic latitudes close to ± 45°. 

The rate of correct detections is lower than the ex- 
pected rate for s t rictly periodic variables as determined by 
lEver fc Mignardl (|2005h because the light-curves of solar-like 
stars are more complex than a simple sinusoidal signal and 
the flux variations due to the rotational modulation can be 
hidden by the intrinsic evolutions of ARs. 

The simulations described in this paper can also be used 
to estimate the rate of correct period detection in other sur- 
veys and, therefore, it could be a very useful tool to estimate 
the compl eteness of rotation period distributions retrieved 
by them. Ilrwin et al.l ([2006) , for example, estimated the 
completeness of the rotation period distributions retrieved 
for the open cluster M34 by running the detection algorithms 
on simulated sinusoidal light curves. The comparison be- 
tween our analysis and that performed bv lEver fe Mignardl 
(|2005h shows that the use of a sinusoidal signal to simulate 
solar-like variability could significantly overestimate the rate 
of correct detections and, therefore, the completeness of ro- 
tation period distribution. 
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